


^ 



> 

o 



X 



INFN-FE 07-94 
July 1994 



Determination of QCD condensates 
^ '• without hadronic spectra 



Murman Margvelashvili^^^'^^Q and Ketino Aladashvili^^^fj] 

(«) INFN sezione di Ferrara, 1-44^00 Ferrara, Italy 



{b) 



High Energy Physics Institute, Tbilisi State University 



cn ; 380086, Tbilisi, Georgia 

O 

A ■ Abstract 

(— I I The bounds on the values of gluon, four-quark and quark-gluon conden- 

sates are derived from the requirement of consistency of the sum rules for 
various correlators of the hybrid current o^ = gd^pGp^u. The upper bound 
for the gluon condensate is found to be less then twice its standard value while 

C^ ' for the four-quark condensates the violation of factorization by a factor 3-4 

is allowed. The value of the quark-gluon condensate is given as a function of 
(^G^) and {{qTq)'^) which for the standard values yields mg = 0.63 GeV^ . 
Our procedure of simultaneous solution of a number of sum rules allows the 
quantitative error analysis and is appropriate for condensate determination 
also in other cases. 
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Introduction 

QCD vacuum condensates, the vacuum expectation values of local operators, in- 
troduced in ref.[0] as essential ingredients of the QCD sum rules (SR) method, are 
the fundamental parameters of theory. They parametrize the properties of QCD 
vacuum and therefore have relation to almost all aspects of the hadron physics. 

Since the introduction of QCD condensates, considerable effort has been de- 
voted to determination of their values, which however up to now has given a little 
effect. The existing theoretical schemes provide only qualitative estimates, while 
semiphenomenological analyses on the basis of QCD SR approach have resulted in 
unacceptably wide ranges of possible values, even for the basic condensates of lowest 
dimensions. 

The most commonly used values of QCD condensates are the so called " standard 
values" (SV) of ref.[|l|. For the gluon condensate this is 

^Gl^ =(330MeV)^ (1) 

TT '^ / St 

as derived from charmonium phenomenology, while for the four-quark condensate - 
the "factorized" value 



(«r«)=>,... = ¥rn<*>' P' 



whith the quark condensate given by < y/aqq >= (—240 MeV)^, and A^(r) - a 
normalization factor depending on the matrix F. 

The mixed quark gluon condensate has been estimated from the baryon sum 
rules to be 






( igqa^uG" —q ) = m^ < qq > (3) 



with ml = 0.8 GeV^ 

In a great number of articles the different types of sum rules (Borel SR 0], 
FESR[|1, Gauss SR [Q etc.) have been applied for extraction of the actual conden- 
sate values from experimental data [^ ^ (here we quote only a few) . Most of them 
use the measured hadronic spectra from e+e~ annihilation or r decay processes. 
These studies however have not led to a consistent picture. The resulting values of 
condensates considerably differ from each other, most of them exceed the SV sev- 
eral times. For the gluon condensate the corresponding factor varies in the range 
(1-6) while for the four-quark condensate it is up to an order of magnitude. As for 
the parameter ttIq, the estimates ranging from 0.2 to 1.1 GeV^ can be found in the 
literature. 

On the other hand there are still arguments against such strong deviations of 
the vacuum condensates from their SVs. So in ref.0 the restrictions on the value 
of gluon condensate have been obtained from the requirement of consistency of 
different SR for the correlators of a hybrid current a^ = gd'~fpGpfj_u. The two point 
functions < 0|T(a^, a^)\0 > and < 0|T('U7^75(i, a^\0 > have been considered and the 
gluon condensate has been bounded by the requirement that these different sum 
rules yield the same result for the pion matrix element < 0|a^|7r(p) >= —ifn^'^PiJ.- 
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Using the factorization assumption for the four-quark condensate, the authors have 
concluded that < -Gf,„ > is forbidden to exceed its standard value by more then 
some 40%. The value of the quark-gluon condensate has been calculated in |^ by 
means the sum rules for the correlator of a^ with the pseudoscalar current. The 
authors, using the SV for the four-quark and the gluon condensates together with 
(5^ = 0.2 from ref. 0, have arrived at the value of eq.(3). 

The mentioned uncertainties of the condensate values lead to considerable ambi- 
guities in the predictions of QCD SR as well as in various other applications. They 
also arise doubts in the consistency of description of the vacuum properties by means 
of local condensates. However considering the situation with condensate determina- 
tion one has to remember that QCD SR is a semiquantitative method and without 
well defined error analysis the results have quite a limited value. Unfotunately not 
all of the attempts of condensate determination meet this requirments. The intrin- 
sic uncertainties of the method related to the truncation of power series and to the 
imperfection of the model spectra are not even parametrized in a fully satisfactory 
way. Hence different types of SR, influenced by these uncertainties in different ways, 
produce the results which seem to be hardly compatible with each other. Presently 
the systematic study of different sum rules with clear statements about related er- 
rors is highly desirable; This will allow either to constrain the values of the basic 
condensates to some reasonable ranges or to show the limits of consistency of the 
whole approach. 

In the present paper we perform the simultaneous analysis of the sum rules 
for different correlators of the current a^ in order to impose bounds on the values 
of the gluon quark-gluon and the four-quark condensates. We are motivated by 
several reasons. First of all let us note that in ref. 0] the gluon condensate has been 
fixed by assuming the SV for the four-quark condensates. On the other hand we 
have quite a reliable estimate based on the r decay data and low energy theorems of 
current algebra, showing that for the definite four-quark condensate the factorization 
approximation is violated by a factor of (2.5-3) 0. So it is desirable to find bounds 
on the above parameters without any ad hoc. assumptiontions. It is also interesting 
to find the allowed region for S^ which is an important parameter determining many 
dynamical characteristics of pion. Our motivation is also to present the framework 
which allows the quantitative error analysis and is more adequate for a systematic 
study of a number of different sum rules. 

Our suggestion briefly is the following. Instead of the standard procedure we 
propose to consider the sum rules at different values of the Borel parameter M^ 
as independent equations with estimated errors. Then we find the minimal x^ 
solution of the system of such equations and test its stability against variations of 
different parameters involved. In section 2 we formulate the problem and present the 
equations used for condensate determination. In section 3 we describe the procedure 
for solution of a system of sum rules. Sect. 4 contains the presentation of results, 
discussion and comments. 



Equations 

In order to impose the bounds on the possible values of different condensates and 
5^ we consider the following set of two point functions of the hybrid current a^: 

I J dxe'"'^ < o|T {u-f^G^,d{x), d-fpG^Mo)) |o >= g.^nHq') + g^g.n^(g') 

tjdxe''^^ < 0\T[u^^^,d{x)Jj,G,MO)) |0 >= g.MTil') + Qf^QuIlTil') (4) 
tjdxe"^'' < 0|T (u-f,d{x),d-f,Gp^u{0)) |0 >= q^UP^q^) 

There are several reasons which make this correlators appropriate for condensate 
determination. First of all, their perturbative expansions start at the two loop level, 
so the nonperturbative terms become dominant and hense easier to determine. In 
contrast to the vector and axial-vector current correlators, in the operator expan- 
sion (OPE) of the two point functions (4) the correlation of gluon and four-quark 
condensates values is negative, which allows to bound their absolute values from 
above. The number of unknown parameters in the sum rules is reduced by the fact 
that imaginary parts of this correlators are contributed by the same physical states. 
These are the pion, oi and higher mass states which in the following will be at- 
tributed to the continuum. In the sum rules we use the usual model spectrum with 
IT and ai as narrow resonances and the continuum equal to the perturbative one, 
starting at some point Sq. After standard manipulations one obtains the following 
SRs for Ilg, 112^11"^ and 11^'^ correspondingly: 



1 /as^,\ ^ , 8{V^sqq).^Q 



-G^) G ' ' '^'^'^ 



72 \ TT /st 9 M2 

= ^ + ^ -p(-m^/M^) - ^M^ (l - (1 + xo + f )e--) (5) 



1 /«.„.\ „ 647r(V^sqq)]Q 



6M2 \ 71 /st 27 M4 

= 1^-1^ exp(-.n^/M^) - ^M^ (l - (1 + .o)e-) (6) 

U7, {V^sqq)]^Q 

27 M2 

= ^-^ eM-ml/M^) - ^M^ (l - (1 + xo + f )e--) (7) 



ml {qq) tt^ (v^') ill) ^ 6' {qq) _ c^s (qq) ^ 
M2 9 M4 M2 37r 



(8) 



where f^ = 133 MeV and fa = 170 MeV|l| are the leptonic decay constants of vr 
and ai mesons correspondingly, m^ = 1260 MeV is the Oi mass and in an analogy 
with the pion matrix element we have e^faS^ =< 0\a^\ai >. We have introduced 
the coefficients 



^)/(v^'L Q^{mr)/{mY)^ 



TT / \ n / St ^ ' ^ ' /«c 



measuring the deviations of corresponding condensates from their standard values. 
Note that Q is some average coefficient for different four-quark condensates which 
appear in the operator product expansion (OPE) of the considered correlators. We 
have Xq = jp where sq is the continuum onset in eqs.(5-7). On the orther hand 
in eq.(8) where there is no ai contribution the corresponding continuumthreshold is 
taken to be lower Sq < Sq . We disregard the anomalous dimensions of the operators 
and consider the corresponding matrix elements at /i^ ~ M^. All the terms except 
the unknown condensate contributions are moved to the right hand sides of the 
corresponding equations. 

A few comments are in order before we proceed further. Equations (5) and 
(6) are taken from 0. In eq.(7) we have corrected the sign of the second term 
in the rhs. which was incorrect in ref.|Q]. This change reduces the resulting value 
of quark-gluon condensate considerably. The SR (8) is not very reliable since its 
different sides show different M^ behavior, however we have added it to the system 
just to illustrate some points of our approach. The second term of the Ihs. of (7) is 
obtained by the factorization of dimension 7 quark-gluon condensate, however due 
to numerical smallness of this term the approximation has practically no effect on 
the results. 

Now let us recall the standard sum rules procedure in order to understand better 
our proposal. Note that QCD sum rules are approximate equations which depend 
on the borel parameter M^. The errors in this equations are mainly caused by two 
independent reasons. At low M^ the inaccuracy of equations is mostly caused by 
the errors resulting from the truncation of OPE series and from the uncertainties 
of theoretical parameters. Whereas at high values of borel parameter the errors 
caused by inaccuracies of experimental (or model) spectra become more essential. 
At the intermeidate values of the Borel parameter there is some interplay of both of 
these errors. If there is a region of M^ where both of them are moderate ("fiducial 
region"), then the physical quantity in question is extracted by equating theoretical 
and experimental parts of SR, as well as their derivatives, at some particular point 
within this region. The reliability of the result is then tested by its stability against 
variation of M^. 

The above procedure is sufficient for deriving qualitative results and estimates, 
however it contains considerable amount of arbitrariness related to the choice of 
fiducial region, determination of the central value and erros of the final result. Be- 
sides it is difficult to compare the results obtained in different sum rules by means of 
such a procedure. The well established quantitative results can be obtained by QCD 
sum rules method only in case if these uncertainties will be somehow eliminated. 

Instead of the traditional treatment we suggest to consider the Borel transformed 
SR (5-8) at different M^ as independent approximate equations in which the esti- 



mates of the expected errors and their correlations should be done. Then the system 
of such equations can be solved in the sense of minimal x^? with account of this er- 
rors. This provides a more natural concept of similarity of different parts of the sum 
rules than the standard procedure. Calculating x^ per degree of freeedom we have 
the measure of consistency of different sum rules entering the system. Then we scan 
all the possible values of the unknown parameters 5"^ and 5l and find the region in 
the three dimensional space of parameters G, Q and tjiq for which the system is 
consistent. This gives us the bounds on the possible values of condensates as well 
as correlations among them. 

Of course such a way is related to many subtleties of the error analysis, however 
this is unavoidable if one wants to controll the precision of the results obtaineds 
by the sum rules method. At least moving in this direction we will investigate the 
capability of the method to provide accurate numerical results. 

Condensate determination 

In order to determine the allowed ranges for the gluon, four-quark and quark-gluon 
condensates we perform the following steps: 

1. Construction of the system of equations. Taking some values of the 
parameters 5^ and 5"^ we evaluate the equations (5-8) at four different values of the 
borel paremeter M^ a = 1,.,4 and build up a system of 16 linear equations for 
determination of the parameters Q, G and mg. We take the M^ points in the range 
0.6 -T- 2 GeV^, however they should not be too close to each other in order to avoid 
strong error correlations among the resulting equations (see below). 

2. The estimate of expected errors. To set the scale of expected errors 
in the resulting equations we first check separately the SRs (5-8) and adopt some 
relative errors Wi for each of them at M^ = 1 GeV^. The eqs. (5-7) are easily 
stabilized by appropriate choice of Sq for reasonable values of condensates and S'^. 
Hence, we have no reason to expect that some important contributions are missing. 
On the other hand eq.(8) exhibits different M^ behavior of its right and left hand 
sides and is clearely less reliable. So, we should take ws > w^^qj. As initial guess 
we assume ^5,6,7 = 0.1 and ws = 0.3 i.e. 10% relative error at M^ = 1 GeV^ for 
eqs. (5-7) and 30% of that for (8). We require that the results shold be stable against 
reasonable variations of these numbers. 

The M^ dependence of errors is taken to be common for all the considered SRs 
and is discribed by the error distribution function 



DiM^) = l^ ^f^^<lGeV^ (9) 

^ ' [I if M2 > 1 GeV^ ^ ^ 

Such a choice of D{M'^) can be justified in the following way; It is known that 
at low M^ the higher order power corrections in the OPE series become essential 
and at some M^ the expansion blows up. So, we assume the pole like behaviour at 
M^ < 1 GeV^ with Mg being some effective convergence radius of the series. On the 



other hand for the higher M^ values the contributions of ai and higher states, which 
become the main source of the errors, remain practically constant, so we assume 
that the error distribution to is uniform for M^ > 1 GeV^. Note that the absolute 
values of SRs decrease with M^ so the above of D{M'^) corresponds to increasing 
relative error both for high and low M^. 

Finally, the expected absolute errors for the equations of our system are calcu- 
lated as 

A,(M2) = uj,D{Ml)A{Ml) (10) 

where Aj(M^) denotes the rhs.s of (5-8). 

3. Solution of the system of equations. To find the allowed regions of the 
condensate values we take different values of the parameters 6'^ and 6^ and for every 
combination we find the solution for Q,G and itlq at which the system is maximally 
consistent. Such a way is dictated by simplicity reasons, since the system is linear 
in condensates. As a measure of consistency of the system we use the minimal x^ 
criterium the use of which we will now try to advocate. 

It is clear, that the errors of our equations are not statistical and should be 
somehow correlated with each other. However these correlations depend on many 
different factors, like unknown higher order terms in OPE, perturbative corrections, 
discrepancies between model and real spectra, possible non OPE contributions etc. 
The knowledge of all these factors would be almost equivalent to the solution of QCD 
for our case, and we are clearly far from this. The best that we can do at present, is 
to consider the errors of different SRs as uncorrelated; As for errors within each of 
the sum rules (5-8), we have to keep the M^ points sufficiently separated to ensure 
that this correlations are small. The effect of our assumption will be less important 
for a greater number of equations, and can be partly compensated by increasing the 
error bars. In other words we suggest to treat the uncertainties caused by a great 
number of unknown parameters as some statistical errors. So, assuming that in each 
of 16 equations the errors are independent we sum them quadratically to calculate 
X^ per degree of freedom. 

We scan all the possible values of the parameters 6'^ and 6^ and for each pair 
find the corresponding values of G, Q and itlq with minimal Xdf- If ^^e resulting 
Xmin < 1 then the system is considered to be consistent and the values of the above 
parameters - allowed. Thus we obtain the allowed region in the three dimensional 
space of parameters G, Q and ttiq. The corresponding region in the plain of (G, Q) 
is shown in fig.l. The inner curve - a) corresponds to Xdfmin — ^ However to be less 
stringent we allow also solutions with Xdfmin — 2 (bounded by the curve 6))which 
corresponds to ~ 15% relative error at M^ = 1 GeV^ in eqs.(5-7) . 

One should keep in mind that for each 5^ and 6^ the parameters G, Q and ttIq, 
giving minimal Xdf solution, have correlated errors which have also to be taken into 
account. For a few solutions this is illustrated by ellipses in fig. 2 . These ellipses, 
for every solution with Xminy include the points for which x"^ ^ Xmin + 1- The area 
covered by all these ellipses is similar to the area with Xmin — 2 bounded by the 
curve b) of fig.l. 
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4. Stability check and parameter adjustment. The resuhs should be stable 
against the variations of the parameters involved in the derivation. So we check the 
stability of the boundaries and individual solutions, and adjust the parametersl in 
order to obtain the maximal stability of the allowed region. 

If the error distribution function is realistic, then the results should not depend 
considerably on the variation of the endpoints Mf and M^. The error distribution 
function (9) with Mq = 0.45 GeV^ just satisfies this requirement. The boundary 
of the allowed region remains practically unchanged when M^ and M^ vary in the 
intervals (0.6 — 0.9) GeV^ and (1.6 — 2) GeV^ respectively. The region shown in fig.l 
is obtained with M^ values taken at 0.7,1.1, 1.4,1.9 GeV^. For comparison in fig. 3 
we have also plotted dependence of the allowed (G, Q) region on the lower endpoint 
Mf for M2 = 0.3 GeVl 

The final result is also insensitive to individual variations of Wi within about a 
factor of 2. However, this is not so for ws and the allowed region feels the erroneous 
eq.(8) unless we take Ws > 0.5. Our final result is not strongly affected by the 
considerable variations of Sq and Sq as can be seen from fig.4. This however is 
hardly surprising in view of the two loop suppression of perturbative contributions. 

Results and discussion 

Thus, we have found the allowed regions for the gluon, four-quark and quark-gluon 
condensates as well as the parameters 6'^ and 6\. We conclude that the system of 
sum rules (5-8) is consistent with the standard values of condensates. From fig.s 1,2 
one can see that the gluon condensate can differ from its SV at most by a factor of 
2 and this happens if the four-quark condensates are close to their factorized values. 
On the other hand the latter are less tied to their SV and can exceed them 3-4 
times. The upper bounds on the values of condensates are more strict than that 
derived from vector and axial channels. This is due to the fact that in the considered 
equations both the gluon and the four-quark condensates have the same sign and 
thus cannot compensate the growth of each other. Note however that we did not 
distinguish different types of the four-quark condensates. 

Since for every S and S^ we obtain a triplet of parameters G, Q and ml, we 
thus have some functional dependence among their possible values, described by 
two dimensional surface in the three-dimentional space of this parameters. This 
dependence is almost linear and can be expressed by the formula: 

ml = (0.246" + 0.3g + 0.09) GeV^ 

which for the SV of condensates (i.e. G=Q=1) yields: 

mo = 0.63 GeV^ 

in difference with eq.(3) 

The allowed regions of condensates correspond to the following ranges of the 
pion and ai matrix elements: 0.1 MeV^ < 5'^ < 0.3 MeV^ and < 5^ < 0.2 MeV^ 
which are thus also bounded by the requirement of consistency of the system of SR. 



The proposed procedure of analysis can be easily extended to other systems of 
two point functions and to a greater number of equations, however in our opinion 
it is already useful for analysis of a single SR, since it allows clear and testable 
assumptions about errors and their distribution. We hope, that this will be helpful 
for further restriction of condensate values or establishing the limits of applicability 
of the QCD SR method. 
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Figure captions 



Fig.l The allowed region in the plane (G, Q) calculated for Mf = 0.7, 1.1, 1.4, 1.9 GeV^ 
The curve a) bounds the region with Xmin — 1 assuming 10% rel. error in eqs.(5-7) 
at M^ = 1 GeV^. b) bounds the region with xLin ^ 2. 



Fig. 2 The allowed region of {G, Q) plane taking into account the errors of the 
solutions with xLin — 1 (inside curve a) of fig.l). The ellipses indicate the areas 
with x^ ^ Xmin + 1 fo^ three different solutions, one of them corresponding to the 

SVs. 



Fig. 3 Allowed (G, Q) regions with xLin — 1 calculated using the error distribu- 
tion function (9) with M^ = 0.3 GeV^; a)for Mf = 0.6 GeV^ b)for Mf = 0.8 GeV^ 



Fig. 4 Allowed regions with Xmin — ^ ^'^^ different choice of Sq and Sg. 
a) - So = 2 GeV=^ s[, = 1.5 GeV^ b) - sq = 3 GeV^ s[, = 2.5 GeV^ 



